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Distortion of the HBT images by meson clouds 

Koichi Hattori*) and T. Matsui**) 

Institute of Physics, University of Tokyo 
Komaba, Tokyo 153-8902, Japan 

We study the effects of mesonic final state interactions on the Hanbury Brown and 
Twiss (HBT) intensity interferometry for mesons in ultra-relativistic heavy ion colhsions. 
Modification of the one-body amplitude of emitted mesons while going through a cloud of 
other mesons is estimated in the semiclassical approximation with a mesonic optical potential 
which incorporates both coherent forward scattering with other mesons and the absorption 
due to the incoherent scattering in the meson clouds. We show how these effects results in 
the distortion of the HBT images. 

Hanbury Brown and Twiss intensity interferometry, originally developed in radio 
astronomy to measure stellar radii/-* and later played a seminal role in the devel- 
opment of quantum optics, has been applied to measure the source size of the 
secondary particles emitted in high energy heavy ion collisions. ^^"^^ Recent analysis 
of the RHIC experiments has shown that the source shape determined from the data 
has a systematic deviation from the predi ction of the hydrodynamical models:^-*"^^ 
this has been called "RHIC HBT puzzle". 



The HBT formula for the correlation function usually used in data analysis 
assumes: ^^-^ 1) random initial phases (incoherent source), 2) factorization of two- 
point source function, 3) neglection of all the interaction between two detected pions 
and the rest of the system after the emission. Each of these underlying assumptions 
may need to be checked in order to find a resolution of the " HBT puzzle" . In this 
work, we investigate the effects of the final state interactions on the correlation 
function, assuming the first two assumptions are valid. 

We focus on the effect of one-body interaction, namely the interacion between 
each of the observed pairs and rest of the system via mean field potential. The 
mutual interaction between the two detected pions is dominated at large separation, 
corresponding to small g, by long range Coulomb interaction, rather than strong 
interaction, and such effect has been incorporated by the well-known Gamow factor 
in the correlation function.^) In the usual hydrodynamical modeling, the particle 
momenta are considered to be frozen on the kinetic freeze-out surface. In the kinetic 
theoretical language, this implies that the collisions between two particles which 
maintains the system in local equilibrium becomes suddenly ineffective. However, 
there still remains interaction between each of the pair and the evaporating particles 
in the vicinity of the emission points. We introduce mesonic optical potential to 
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describe such effects: tlie real part of the potential describes a coherent forward 
scattering of the particle with the other mesons, while the imaginary absorptive part 
incorporates the effect of incoherent scattering with other individual mesons. We 
examine how extracted HBT images are distorted due to the modification of the 
one-body amplitudes in the mesonic optical potential. 

In the standard picture of ultrarelativistic nuclear collisions, the pion source 
extends along the direction of the motion of colliding nuclei in an approximately 
boost invariant fashion. We follow this picture, together with cylindrical symmetry 
of the collision volume, and concentrate on the quantum evolution of a group of 
particles appearing in a small central rapidity bin —\Ay < y < on the two- 
dimensional transverse plane in their center-of-mass frame. 

Two particle momentum correlation function of identical particle which we are 
concerned is given by 

where ^2(^1; ^2) is the joint probability of detecting a pair of the same kind of pions 
with momenta ki and k2, and Pi{k) is the probability for a detection of single pion. 
They may be given in terms of the asymptotic form of the matrix elements of the 
density matrix: Pi{k) = {k\p{oo)\k) , ^2(^1,^2) = {kik2\p{oo)\kik2). Here we 
adopt the Schrodinger picture in which the time evolution of the density matrix is 
given by p{t) = U (t) p{0)U^ {t) so that we may write 



Pi(fc) = lim / dxidx2{k\U{t)\xi){xi\p{0)\x2){x2\U^{t)\k) , 

t-^ooj 

= lim / dxidx2(pk{xi,t)po{xi,X2)(pl{x2,t) (2) 

t^oo J 

-^2(^1,^2) = lim dxidx2dxidx2{kik2\U{t)\xix'i) 

t— »oo J 

x{xi,X2\p{0)\x[,x'2){X2x'2\U^{t)\kik2) 

= lim / dxidx2dx[dx'2^ki,k2i^i^^-2'-'''^)Poixi,X2;x[,x'2)^l k2i^[,x'2,t) 

t— >oo J ' 

(3) 

where ipk{x,t) = {k\U{t)\x) is the amplitude that a particle emitted at x is detected 
with momentum k at time t, while !f'fci,fc2(^i5 ^2, = {kik2\U{t)\xiX2) is the two- 
particle amplitude that two particles are emitted at xi and X2 and detected with 
momenta ki and ^2 at time t simultaneously. We have introduced the one-particle 
reduced density matrix, pi{xi,X2) = (a;i|p(0)|a;2)5 and the two particle reduced den- 
sity matrix, ^2(2=1 ) = {xi,X2\p{0)\x'^,X2) which obey, due to the bosonic 
symmetry of the two particle states {\xi,X2) = |a;2,£Ci)) , 

P2{xi,X2;x[,X2) = P2{X2,Xi;x[,x'2) = p2{Xi,X2;X2,x[) (4) 

The basic objective of the HBT interferometry is to extract information of the density 
matrix from the observed correlation function C{ki,k2)- 
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If the particles do not interact after t = 0, the time evolution operator turns into 
a trivial phase factor and the single particle amplitude becomes just a plane wave: 



Ht\ 



-iE^t 



{k\x) = e-*('=-'^+^fc*). This results in 



(5) 



where f(x,k) with x = {xi + X2)/2 is the one-body phase space distribution func- 
tion or the Wigner function at t = 0. Similarly, in the absence of the final state 
interaction, the symmetrized two particle amplitude becomes, {kik2\U{t)\xiX2) = 
^g-i(fci.a;i+fc2-a;2) +e-i{fci-a;2+fc2-3:i)] ^-i{Er+E2)t _ j£ -^^^^^^ ^-^is expression into 

Eq ([3]) and further assume the factorization of the two particle density matrix, 
Po{x-i,X2]x{,x'2) = \ [po{xi,x\)po{x2,x'2) + Po{x2,x'^)pq{xi,x'2)], along with the 
symmetry Q, we find 
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P2{ki,k2) = Pi{ki)Pi{k2) + 



dxf{x, k)e 



iqx 



(6) 



with k = ^{ki + k2) and q = fci — fc2- For a factorized form of the phase space 
distribution f{x,k) = p[x)Pi[k) with J" dxp{x) = 1, this gives the familiar form of 
the correlation function, 
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C{q,k) = l + r]{ki,k2) 



dxp{x)e 



tqx 



(7) 



with a correction factor r][ki,k2) = (fc)/Pi(fci)Pi(/c2), which exhibits limiting 
behaviors: C(0, k) = 2 and C{oo, k) = 1. This formula has been used to reconstruct 
the source distribution function from the measured two particle momentum correla- 
tion. We note that the lesson we learn from this simple exercise is that the phase 
interference of the symmetrized two-particle amplitude plays essential role to create 
the two particle momentum correlation. 

We now consider how this result is modified in the presence of the final state 
interaction. We shall study the change of the single particle amplitude {k\U{t)\x) 
as well as the two particle amplitude {kik2\U{t)\xiX2) by the action of the time 
evolution operator U{t) = e~*^*. Here we replace the many-body Hamiltonian H 
by one-body Hamiltonian for a non-relativistic particle propagating in an one-body 
potential. Without mutual two-body interaction of pion pair, two-body amplitude 
in ([3]) is expressed by the symmetrized products of one-body amplitudes : 



{kik2\Uit)\xiX2) 



i {{ki\U{t)\xi){k2\U{t)\x2) + {k,\U{t)\x2){k2\U{t)\xi)} (8) 

We therefore only need to compute the distortion of the single particle amplitude. 

We recall that ipkix, t) = {k\U{t)\x) denotes the amplitude for a particle emitted 
at X at time i = being detected with asymptotic momentum k at sufficiently large 
time t. Formally, it may be interpreted as a wave function which obeys the static 
one-body Schordinger equation on the transverse plane: 



1 

2m 



+ V{x) 



iPk{x,t) = Ek^Pk{x,t) 



(9) 
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with Ek = k'^/2m. If we assume cylindrical symmetry for the collision volume, this 
Schrodinger equation in two space dimension may be reduced to a one-dimensional 
one in a "central" potential y(|a;|) and may be solved by the standard WKB semi- 
classical approximation. 

Here we adopt more intuitive approach to the problem motivated by the station- 
ary phase approximation of the path integral expression of the amplitude, exploiting 
an analogy to geometrical optics. We write 

{k\U{T)\x) = J dx'e-'^-"^' {x'\U{T)\x) (10) 

and introduce an Ansatz 

(a;'|;7(r)|a;) « A(a;', a;; r)e^^=(^''^;^'°) (11) 
for the propagator. We work with the classical action integral 

Sc{x2,t2;xi,ti) = / dtL(x,x) (12) 
Jti 

along the classical trajectory {x{t),p{t)) specified by the initial position Xi of the 
particle at time t = ti and the final position X2 at time t = t2. The initial and 
final momenta of the particle are given by pi = —'S/iSc{x2,t2'-,xi,ti) and p2 = 
^ 2Sc{x2,t2]Xi,ti) with the nablas Vi^2 operating on xi^2^ respectively. For free 
particle motion 

5f-(^',r;^,0) = ^^^i^^^ (13) 

so that the integration over x' , assuming that the prefactor A{x^ x' ^ T) is a constant, 
evidently give the previous resuhs: {k\U{T)\x) = e-**'-^-»-^fc'^. This motivates us to 
write: 

{k\U{T)\x) = A{x, fc)e^w^(^.fc)-i-E^fe^ = ip^{x)e-'^''^ (14) 

where we separated the phase factor into an explicitly time-dependent piece E^T 
and T-independent part W{x, k). Evidently, for free particle, W{x, k) = Wo{x, k) = 
—k • X. 

To construct the phase factor W{x, k) in the presence of interaction, we decom- 
pose the classical action into two parts: 

Sc{x', T; X, 0) = Scix', T; x" , T') + Sc{x" , T'; x, 0) 

where {x" ,T') is chosen at any point x"(T') on the classical trajectory. If we choose 
a point sufficient far away from the interacting region, we may set 



Sc{x',T;x",T') = So{x',T;x",T' 



mix' - x"Y 



2{T - V) 

Putting this into (jlip and performing the integral over aj', we obtain, 

{k\U{T)\x) = A{x',x, T)e-^^^(^-^')-*'^-"+^^=(-"'^''-'0) (15) 
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Comparison with (I14p gives the formula: 

W{x, k) = EkT' -k-x" + Sc{x",T'; x, 0) (16) 

and A{x,k) = A{x' ,x;T). We emphasize that {x",T') is just a dummy coordinate 
which can be arbitrarily chosen on the classical trajectory as long as T' is taken to 
be sufficiently large so that is a function only of x and k. 

This procedure determines W{x,k) uniquely if the classical trajectory corre- 
sponding to {x, k) is given. There are cases, however, that there are more than one 
classical trajectories which correspond to the same x and k, similar to appearance 
of caustics in geometrical optics although time is reserved in our case. We will see 
this actually occurs in the presence of interaction. 

In the presence of the imaginary part V2{x) of the optical potential the prefactor 
may be given 



A{x, k) = a{k) exp 



dtV2{x{t)) 



(17) 



where the integration is along the classical trajectory determined by the real part of 
the optical potential Vi{x) for given x and k. This integral is negative for absorptive 
potential V2{x) < 0. 

We first consider how Pi{k) is modified by the interaction. Using the center- 
of-mass coordinate x = ^{xi + X2) and relative coordinate r = xi — X2, we have 
^kixi) = ^Pk{x + ^r) ~ exp [i^r • VW{x,k)] ipk{x), and ipk{x2) = ipk{x - ^r) ~ 
exp [— i^T" • VW{x, fe)] ipk{x), ignoring the higher order derivatives of W{x, k) with 
respect to x, as well as the derivatives of the prefactor A{x, k) with respect to x. 
This procedure may be justified for our choice of Vi{x) which is a smooth function 
in the source region. Inserting the above results in ([2]) we find 

Pi{k) = j dxdr j ^|^e-^^-('=-^^("''=))/(a^,fc)^'(a;,fc) 

= a^{k) j da;/(a;,p(a;,fc))e-2^(^'*') (18) 

where p{x, k) = ^W{x, k) is the initial momentum of the particle when it is emitted 
and 7(2;, k) = — dtV2 measures the amount of absorption. 

This result is precisely what is expected by our intuition. The real part of the 
optical potential causes the shift of the momentum of the observed particle from the 
original momentum at the time of emission by acceleration by the potential field. 
This is purely classical effect which may be obtained from the classical treatment. ^"^^ 
In the quantum description, the information of the momentum shift of the particle 
is encoded in the phase shift of the single particle amplitude in proportion to x. The 
flux attenuation factor e~^'^(^'*') is due to the absorption (including elastic scattering 
of the particle) on the way out of the emission point which may also included on 
purely classical consideration. 

We now study how this phase shift results in the change of the two-particle 
momentum correlation through the change of interference which arises from the 
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symmetrization of the two-particle amplitude. The two particle joint probability ([3]) 
may be evaluated with the same technique. Inserting 

and similar form for l^'^^ (x'^, X2,t) in ([3]), and then assuming again the factorization 
of two-body reduced density matrix, we obtain 

P2{ki,k2) = Pl{k^)Pl{k2) + \Fiki,k2)f (19) 

where we have introduced the integral 

F{ki,k2) = j dxidx2Pi{xi,X2)(pkAxi)fk2i^'i)^ (20) 
to express the interference term. Using the Ansatz (jlip and ([1] 



F{ki,k2) = a{ki)a{k2) J dxidx2p{xi,X2)e 



i{W{xi,k2)-W{x2,ki))g-'r{xi,k2)-'r{x2,ki) 



(21) 



We approximate the exponents by W{xi, ki) — W{x2, ^2) = W{x + ^r, k + |g) — 
W{x — ir, k — ig) =~ r ■ VW{x, k) + q ■ VkW{x, k) and 7(2:1, ^2) + 7(2^2, ^1) — 
2'j{x,k) , retaining only the leading terms in the expansion in q since we expect 
to find a significant momentum correlation only for small q less than the inverse of 
the linear dimension of the source. Inserting p{xi,X2) = (2~p / dpf{x,p)e~'^P''^ and 
then performing integral over r, we find 

F{ki,k2) = a{ki)a{k2) J dxf{x,p{x,k))e-"'-^''^^''^''^e~'^^^''^''^ (22) 

where p{x, k) = VWlx, k) is the initial momentum of the particle. 

For free particles, we may set W{x, k) = Wo{x, k) = —x k, a{k) = 1, 27(03, k) = 
and k = p so that (f22l) is reduced to 

F{ki,k2) = J dxf{x,k)e''^-'^ (23) 

which reproduces ([6]). 

To see how the interaction modifies the source image we write 

W{x, k) = Wo{x, k) + 6W{x, k) = -x ■ k + 6W{x, k) (24) 

where 5W{x, k) is the phase shift by the mean field interaction. Since VfcVF(a;, k) = 
-X VubWix, k), this gives ([22]) as 

F{ki,k2) = a{ki)a{k2) J da;/(a;,p(a;, fc))e^''-(^-^'='^^(^''=»e-2''('^''=) (25) 

Note that the phase shift 6W{x, k) caused by the real part of the mean field potential 
resulted in the apparent shift of the emission point by 5x = —Vk^W^x, k) while the 
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we can write 



imaginary part effectively cut off tlie contribution of tlie source from deep interior 
and the back side. Transforming the integration variable x to 

x' = x-VkSW{x,k) (26) 

F{ki,k2) = a{ki)a{k2) I dx'f,s{x', k)e'^-'^' (27) 
with 

/eff(a;', k) = J{x, x'; k)f{x,p{x, fc))e-2^("''=) (28) 

where J{x,x';k) = d{x,k)/d{x' ,k) = [d{x' ,k)/d{x,k)]~^ is the Jacobian of the 
coordinate transformation {x,k) — > {x',k) defined by (j26p and x on the right-hand 
side is understood as a function of x' and k. The correlation function is now given 
by 



C(fci,fc2) = 1 + 



dx'p^s{x',k)e'''-''' 



2 



(29) 



with the effective source distribution defined by 

~ ( ' i,\ — fcsjx', k) 

Note that the momentum dependent prefactors 0(^1)0(^2) cancel out by the nor- 
malization. The absorption factor e~'^"'^'^'^^ gives a weight on the distribution of the 
emission points toward the side of the observation of the particle. 

For numerical computation we have chosen here a schematic pion optical po- 
tential with two-range Gaussian shape for the real part, Vr{x) = Vie~^ /^•^i + 
V2e~^l^l~^^ /2-'*2j and one-range Gaussian for the imaginary part, Vi{x) = —Vie~^ Z^-**!^ 
where the shorter range Ai of the potentials is taken to be the same as the parameter 
for the Gaussian source distribution, p{x) = po^"^ while A2 gives the thickness 
of the meson hallow surrounding the dense meson source. We have chosen Ai = 5 
fm, A2 = 5 fm, C = 10 fm, Vi = 10 MeV, V2 = 2 MeV, Vi = 0.1 MeV. We also used 
thermal Bose-Einstein distribution of the initial pion momentum with temperature 
To = 140MeV. 

The classical trajectories are computed numerically with the initial condition 
{x,p) until the asymptotic momentum k is reached. Since the total energy is con- 
served the magnitude of k is determined by the initial value of \p\ and Vi{x). The 
integration gives the angle of deflection with respect to the initial direction of emis- 
sion p. Along each trajectory the action integral is computed numerically; this 
determines the phase factor W{x, k) as a function of {x, k). Computing many such 
trajectories with different initial conditions, we have a good sample of trajectories 
and the phase factor W{x,k). We use these data set to compute the derivatives 
of 5W{x,k) with respect to \k\ and 9k. by interpolation which are related to the 
apparent location of the emission points in Cartesian coordinate, x' = {x',y'), where 
x' = X — is the apparent location of the emission point in the direction of the 
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Fig. 1. Contour plots of the effective source distribution p^a.ix' ,k) defined by (31) at fc = 
lOOMeV/c. Left (right) columns correspond to results with repulsive (attractive) potential. 
Lower panels include effect of absorption. Vertical (horizontal) coordinate denotes the shifted 
outward (sideward) location of the emission points whose plot range is distorted by the non- 
linear coordinate transformation. 
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Fig. 2. Contour plots of the effective source distribution at fc = 150MeV/c. Each panel corresponds 
to the same choice of the potential at the same location in Fig.l 



final meson momentum k (outward direction) , while y' = y — \ is the apparent 
location perpendicular to k (sideward direction) where 9^ is the angle of k. 

We show in Fig. 1 contour plots of the effective source image at A; = lOOMeV/c. 
The upper two panels give the results with no absorption (Vi = 0) while the lower 
panels are computed with Vj = 0.1 MeV. The vertical and horizontal coordinates 
gives the apparent location of the emission points in the outward direction (direction 
of k) and the sideward direction (direction perpendicular to k) respectively. We first 
compare the case for repulsive interaction Vi = 10 MeV with no absorption (upper 
left panel) and the attractive interaction with no absorption (upper right panel). We 
found that the real part of potential affects on the extension of source in both outward 
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and sideward directions. The repulsive interaction leads to elongation of the effective 
source image in the outward direction, while the attraction tends to shrink the source 
in this direction. On the other hand, the repulsive force leads to shrinking sideward 
source extension, while the attractive force gives the sideward extension stretched. 
These effects are seen as a geometrical effect of either stretching or compressing the 
scale of each coordinate. The distortion of the original image is weakened as the 
momenta of the two particles increases as seen in Fig. 2 at A; = 150MeV/c. 

Our results qualitatively agree with the results of two other groups. Pratt 
and Cramer et al. suggested that this change of the apparent source size in the 
sideward direction may be interpreted as due to the refraction or the lensing effect in 
geometrical optics. In our analyses, however, the apparent shift of the emission point 
is caused by the k dependence of the phase shift 5W{x,k), through Eq. (26): the 
apparent shift of the emission point dx = —VkSWlx, k) arises due to the difference 
of the action of the potential for particles traversing along two adjacent trajectories 
starting at the same point x ending up with different final momenta. 

To make this issue more quantitative, we examine the phase shift using Glauber- 
type approximation which assumes the straight-line trajectories in the interaction 
regions so that the phase shift is given by the simple formula which may be written 
for F << fc^/2m, 

^ fu{T) 

5W^Glauber(a;,fc)~-— / V{^A^^+P)du (31) 

where b = \J x'^ — {x ■ k)'^ jU^ is the "impact parameter" of the trajectory, and ii(t) = 
xit) ■ k/\k\ is the position of the particle along the trajectory at time t. This phase 
shift depends on 9k since b is related to fc by 6 = ^/x{T)^ - {x{T) ■ kf/k'^ at large 
time T. The change of the HBT image calculated by this formula is shown in 
Fig. 3 for the same cases as in Fig.l. It is seen that the "Glauber approximation" 
qualitatively reproduces the same results. This implies that the deflection of the 
classical trajectory in the source region is not essential for the distortion of the 
HBT images; rather, the change of the relative momentum of two particles via the 
difference of their phase shifts is the origin of the distortion of the source image. 

The inclusion of the absorption, however, diminishes all these interesting effects 
of the mean field interaction, by effectively cutting off the contributions from deep 
interior of the source as well as the other side of the source. The result is somewhat 
similar to the effect caused by the attractive interaction, namely the source image is 
effectively stretched in the sideward direction. 

Although we still need to make improvements in our calculations (relativistic 
treatment, time-dependent source structure, ^^^'^^^ more realistic pion optical poten- 
tial, etc.) before confronting the experimental data, these mean field effects may 
play some role in reducing the discrepancy between the data and hydrodynamic 
simulations at small values of k. 

In conclusion, we have studied the effect of the final state interaction in the 
meson clouds on the HBT interferometry in heavy-ion collisions and have shown 
that the final state interaction causes a significant distortion of the source images at 
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y' :sideward(fm) 




Fig. 3. The effective source distributions ai k — f OOMeV/c in tiie Glauber approximation computed 
witli tfie same potentials as in Fig.l 



small k through the change of the phase shift in the single particle amplitude and 
the absorption effect. More detail account of this work will be reported elsewhere. 

Acknowledgements 

We thank Gordon Baym, Tetsufumi Hirano and Koichi Yazaki for helpful con- 
versations on the related works, and Hirotsugu Fujii for calling our attention to the 
reference. This work is supported in part by the Grants-in-Aid of MEXT, Japan, 
No. 19540269, and Global COE Program "the Physical Sciences Frontier", MEXT, 
Japan. 



R. Hanbury Brown and R. Q. Twiss, Nature 177, 27 (f956) 
R. J. Glauber, Phys. Rev. Letts. 10, 84 (f963) 

F. B. Yano and S. Koonin, Phys. Letts. B78, 556 (1978) 

M. Gyulassy, S. K. Kauffmann, L. W. Wilson, Phys. Rev. C20, 2267 (1979) 

G. Baym, Acta. Phys. Polon. B 29 f839 (1998) 

C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 87 082301 (2001) 

K. Adcox et al. (PHENIX Collaboration), Phys. Rev. Lett. 88 f92302 (2002) 
T. Hirano, K. Tsuda, Phys. Rev. C 66 054905 (2002) 

M. A. Lisa, S. Pratt, R. Sohz and U. Wiedemann, Ann. Rev. Nucl. Part. Sci. 55 357 (2005) 

S. Adler et al. (PHENIX collaboration), Phys. Rev. Lett. 98, 132301 (2007); S. Afanasiev 

et al. (PHENIX collaboration), Phys. Rev. Lett. 100, 232301 (2008) 

P. Danielewicz, S. Pratt, Phys. Rev. C 75 034907 (2007) 

M. C. Chu, S. Gardner, T. Matsui, R. Seki, Phys. Rev. C 50 3079 (1994) 

J. D. Bjorken, Phys. Rev. D27, 865 (1983) 

Similar analysis has been done for the long range Coulomb interaction in G. Baym, P. 
Braun-Munzinger, Nucl. Phys. A 610 286c (1996) 

G. Cramer, G. Miller, J. Wu and Jin-Hee Yoon, Phys. Rev. Lett. 94 102302 (2005) 

S. Pratt, Phys. Rev. C 73 024901 (2006) 

G. Bertsch, G. Brown, Phys. Rev. C 40 1830 (1989) 

D. Rischke, M. Gyulassy, Nucl. Phys. A 608 479 (1996) 



